Restoration of the residue factorizability in the bound-state pole by 
instanton-antiinstanton configurations 
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The instanton-antiinstnton contributions to the qq bound state pole in the four-point Green 
function in the Schwinger Model are calculated. It is shown that these configurations, thanks to the 
cancellation of all unwanted terms, are responsible for the restoration of the perfect factorizability 
of the residue. 
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In gauge field theories the interplay between gauge 
group and space-time dimension may lead to the appear- 
ance of nontrivial topological effects, resulting from ho- 
motopy properties. The most important example of such 
a theory is Quantum Chromodynamics, where specific 
(instanton) gluon-field configurations play a significant 
role in the formation of the intricate vacuum and are 
related to various effects like chiral symmetry breaking 
and nonzero value of the quark condensate, provide solu- 
tion of the 17(1) problem or influence quark interactions 
and meson or baryon correlators, to enumerate only few 
examples 0, 0, [1 0, 0. To clarify these issues much 
work has been devoted to the study of instantonic ef- 
fects in model theories and particularly in the Schwinger 
Model [10,@, BESOM 51, which is highly nontrivial 
and — due to its similarity to QCD in many aspects — 
of real physical importance [l4| • The significance of this 
model for the investigations in hadronic physics should 
be emphasized especially as it allows exact and analytic 
results for a variety of interesting quantities. 

In our previous paper [l5| we concentrated on the 
Bethe-Salpeter wave-function for the qq bound state, 
picked out from the polar term of the two-fermion Green 
function. This Green function was found exactly both 
with and without instanton contributions fl3l . [la . [l7| . 
The incorporation of higher instanton sectors (i.e. for 
instanton number k = ±1,±2) turned out not only to 
modify the form of the Bethe-Salpeter function, which 
acquired additional terms, but also to destroy the factor- 
ization property of the residue in the bound state pole. 
This factorization was still maintained on the level of the 
S matrix, but proved to be spoilt in the Green function. 
The other instantonic effect was the appearance of non- 
polar, branch-point singularities, for P 2 = /i 2 , where /z is 
the invariant mass of the bound state. The goal of this 
Brief Report is to show that the obtained nonfactoriza- 
tion of the bound-state pole residue may be cured, if one 
includes into consideration instanton-antiinstanton {I A) 
configurations. 

The calculation of the I A effects is more challeng- 



ing than ordinary instanton calculus. Such a config- 
uration formally bears k = topological number and 
there are only approximate fermionic zero modes of the 
Dirac operator, contrary to the latter case, where Atiyah- 
Singer index theorem guarantees the existence of true 
zero modes [l8| . However, when the separation between 
instanton and antiinstanton becomes very large, these 
quasi zero modes approach real ones. Happily this is just 
the situation, we are interested in. Picking up the polar 
term out of the coordinate-space two-fermion Green func- 
tion requires Fourier integration over infinite space, and 
the leading contribution is determined by the behaviour 
at space-time infinity. One can then make instanton- 
antiinstanton separation arbitrarily large. 

Following 0, [IH we substitute into the path integral 
in the generating functional 



the following form of the gauge potential 

A»(x) =A^{x)+e tlv d v b{x) , 



(1) 



(2) 



with e^ u being the antisymmetric symbol. The integra- 
tion, which will be led in cuclidean space, is now per- 
formed over A {Q ^ and is restricted to the instanton sector 
k = 0, which requires the topologically trivial behaviour 
of the gauge field at infinity. The substitution Q is a sim- 
ple shift and does not influence the integration measure. 
The background function b is chosen below for our con- 
venience. C is the standard Schwinger Model Lagrangian 

C(x) = *(as) (iY*0„ - g^A^x)) ^{x)-^F^(x)F flv (x) . 

(3) 

with a gauge fixing term, if needed. 1 Since we are in- 
terested in the I A configuration, b will be taken as a 
simple sum of the fields corresponding to pure instanton 
and pure antiinstanton at shifted positions 0, Hfl| ■ The 
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1 Green functions are, naturally, gauge dependent. The following 
formulae will be given in the Landau gauge. 



2 



function b has, therefore, the following (euclidean) form 



b (+-)( x ) 



i_ f{x + R/2) 2 + p 2 

i , f(x-R/2) 2 + p 2 
— In ' 
2ff 



(4) 



where (H — ) refers to the I A configuration. For AI one, 

we will have bi \-)(x) = —bt^ )(x). The separation of I 

and A is set equal to R and we will be interested in the 
limit R — > oo. 

The functional integration over A^^ in is simple. 
First we have to gauge away the term gty A 1 - ^ from 
the Lagrangian by the appropriate redefinition of fermion 
fields (which leads to the appearance of the gauge-boson 
mass term with p? = g 2 /n). Next A^ ^ appearing in 
source terms (because of the above redefinition) is re- 
placed with the functional derivative over external cur- 
rent J p and may be driven out from under the integral. 
The remaining integral becomes then Gaussian and easy 
to be taken. The details of this calculation are given 
in [l3j and are not dependent on the specific choice of 
the function b, hence there is no need to repeat them 
here. 

The next step is to perform fermion integrals. They are 
more complicated but, fortunately, the majority of the 
work has already been done in [l3| . We limit ourselves to 
the four-point Green function, which is the fourth-order 
coefficient of the expansion of (TT]) in powers of fermionic 
sources r\ and rj (i.e. it is G from the term TjrjGrir)). This 
four-point function was found exactly with all the instan- 



tonic corrections (k — 0, ±1, ±2), as well as the form of 
the qq bound state pol e and the appropriate formulae 
may be found in [UGl. What 

we would like to concen- 
trate on below, is the sector k = with the background 
configuration ([4]), when R approaches infinity. This cal- 
culation is somewhat similar to that for k = ±2 in the 
sense that we have again two (quasi) zero modes, which 
(both) have to appear in the Green function due to the 
rules of Grassman integration (cf. formula (51) of [HI]). 
We are only concerned with the leading asymptotic term, 
so these two quasi zero modes may be taken as that of the 
pure instanton centered on —RJ2 and that of the pure 
antiinstanton centered on R/2 0, OH HH • Obviously, for 
finite R they are not exact. 



X2(x) = 



1 



1 



2^ \(x + R/2) 2 +p 2 



1 



1 



2^ \(x-R/2) 2 +p 2 



1/2 



1/2 



(5) 
(6) 



They are eigenfunctions (with eigenvalues ±1) of the chi- 
rality operator 7 s (our conventions as to the 7 matrices 
are given in [HI). 

Now our previous calculation can be applied with only 
obvious modifications and with the additional difference 
that our formulae will become precise only in the limit 
R — > 00. Referring the Reader to [13| and without going 
into details, we can write down the formula for the I A 
contribution to the Green function in question 



1 



Xla(xi)xtc{x3)X2b( x 2)xtd( X 'i)e 



-ig(b {+ _ ) (x 1 )+b ( + _ ) (x 3 )-b ( + _ ) (x 2 )-b ( + ^ ) (x4)) 



A1A2 

+ ^ik^ki^k^K 9 ^ 

ig 2 l0(x 1 -x 3 )+f3(x2-x 4 )-f3(x 1 -x 4 )-f3(x 1 -x 2 )-0{x2-x 3 )-!3{x3-x 4 ,)] _ I C <-> d 



(7) 



x e 



x% <-> X4 



Here a, b, c, d are spinor indices, the function (3 was de- 
fined in [131 ] and Ai^ are eigenvalues for two quasi zero 
modes of the Dirac operator given below in ((5|) and ap- 
pear in denominator due to the normalization of the 
vacuum-vacuum transition amplitude to unity in the 
k = topological sector (other sectors do not contribute 
to this amplitude if fermions are massless, since the tun- 
neling between various topological vacua is suppressed 
by zero eigenvalues). We would like to stress here again 
that Xi 2 are n °t eigenvectors with eigenvalues A1.2, but 



approach them when R tends to infinity. Neither quasi 
zero modes nor their eigenvalues Ai^ have to be known 
exactly, but it is sufficient to know their asymptotic be- 
haviour (naturally Ai^ tend to zero). 

Now we have to estimate these eigenvalues for large R. 
In this limit the true quasi zero modes can be chosen to be 
arbitrarily close to Xi 2- Therefore, following (5L I22I l23l| 
we assume, that the subspace in question is spanned just 
by Xi and % 2 - The euclidean Dirac operator in the I A 
background has the form 



3 



D 



I A 



id 2 - <9i 
id 2 +dx 

I ( 



1 



(a; + R/2Y 


R 1 /2 + i(x 2 







x 1 +R 1 /2-i{x 2 +R 2 /2) 



.2 \x 1 +R 1 /2 + i(x 2 +R 2 /2) 
Xx-R l /2-i(x 2 -R2/2) 

R2/2) 







(8) 



In the vector space with the scalar product defined as 
</ls> = / d 2 xf + (x)g(x), it is a hcrmitian operator with 
only off-diagonal matrix elements 





u 1.\ 



n IA 





(9) 



^YhPW^ where indices 1 



and eigenvalues equal to ± 

and 2 refer to the functions ^ and ([6|). The diagonal 
elements disappear since §E§ inverts the chirality and Xi 2 
are chirality eigenvectors. 

One can easily verify, that D\\ = D f ±, which means 
that the eigenvalues are real, as necessary for the hermi- 
tian operator. These eigenvalues should be substituted 
for A1.2 in frj). but what we need, is only their asymptotic 
form. For instance we have 



D 



I A 



1 

2tt 



= ^ / d 2 



d*xxl{x)D 1A X x{x) 

X\ 



(10) 



ix 2 



(x 2 + p 2 ) 3 / 2 ((x + R) 2 + p 2 ) 1 / 2 



where the integration variable has been shifted by R/2. 
This integral may be performed using the method of 
Feynman parameters, usually applied in the calculation 
of Feynman diagrams, and we obtain 



D 



IA 



1 



1 v /q(l-a)(i?i-zfl 2 ) 
a R 2 a(l -a)+p 2 



Ri + iR 2 
(11) 

which leads to the asymptotic form for the product of 
eigenvalues Ai ■ X 2 ~ -gj-. 

Now we can come back to To find the polar contri- 
bution, we introduce new variables, as it was done in [l5[ 



X = -{xi +x 3 ) , 



x = x x -x 3 , 



Y = -(x 2 +x 4 ), y = x 2 -x 4 - (12) 

The translational invariance should manifest itself 
through the dependence of G on x\ — Xj 's only or, in new 
variables, on Z = Y — X, x and y. This is not visible 
in |[7J) since, for finite R, our formulae are only approx- 
imate ones. Because we are interested in the ^-channel 
singularity, the pole corresponding to the bound state 
should be found in the complex plane of P 2 , where P is 
the two-momentum canonically conjugated to Z . The Z 
dependence of the expression is then crucial. 



Let us begin with the first term in the square brackets 
in 1(7} ■ Its matrical structure, according to ([5]) and ([6]), 
has the form (1 — 7 5 ) ac (l + J 5 )bd- To reveal the de- 
pendence on Z we first have to simplify the expression 
exploiting (|3|) , (JSJ) and ^ , which leads to the explicit can- 
cellation of several factors coming from \ and \i 2 - 

Next the Z-dependent function should be isolated of the 
last exponent (that containing j3 functions). Up to the 
opposite sign, this was actually done in [151 ] (cf. formu- 
lae (26), (27) and following ones). The Z dependence is 
there twofold: firstly in the factor (with altered sign) 



ti 2 Z 2 



exp[ln(-^ 2 Z 2 /4)] 



(13) 



and secondly in e~ lpZ , which (for the polar term) simply 
shifts the Fourier variable P. The rest of the expression 
(i.e. the pole itself, although, before taking the Fourier 
transform, in the variable p ) may be derived from the 
formula (30). 

Gathering all essential factors we see, that the func- 
tion (still written in euclidean space) to be Fourier trans- 
formed is 



R 2 Z 2 



Ox X Oy (Jy 



(14) 



where cr± = {{X ± x/2 - R/2) 2 + p 2 ) 1 ' 2 and af = ((Y± 
2 //2 + J R/2) 2 +p 2 ) 1 / 2 . 

For the limit of infinite R the relative variables x and 
y become unessential and may be omitted. X and Y, in 
turn, can be rewritten as 

X = (X + Y)/2- Z/2 
Y = (X + Y)/2 + Z/2 . 

The polar term is connected only with the infinite Z in- 
tegration, so (X + Y)/2 may again be disregarded as 
R — > 00. This leads to the simplified expression 



\m 2 z 2 



(Z + R) 2 + Ap 2 



(15) 



Now we are in a position to perform the Fourier inte- 
gration over Z , together with the limit over R. They 
should be done in such a way that Z and R are always 
close to each other ((Z — R) 2 « R 2 ), so one can substi- 
tute Z = R + w in (Tj"5")) , take the asymptotic term as R 
goes to infinity, and finally perform the Fourier integral 



4 



over new variable w. The Fourier factor under the inte- 
gral is e t ^ F ~ p ' z . Analyzing (fT"5)) one sees that asymptot- 
ically the integrand function becomes unity and, there- 
fore, the leading term of the Fourier transform is simply 
(2-k) 2 5^(P-p). 

In an analogous way one can easily verify that the sec- 
ond term in the square brackets in ([7]) does not con- 
tribute to the pole. The term obtained by the anti- 
symmetrization, in turn, contributes only to the pole in 
the u-channcl. 

We can now gather the whole expression for the I A 
configuration 



iG£XP;x,y)=- 



fi 2 cos(Px/2)cos(Py/2) 



8tt 



P 2 - /i 2 



(16) 



x e 2lE e l92Wx)+my) \l - 7 5 ) <g> (1 + 7 5 ) , 

where b.s. stands for 'bound state'. The coefficient /k 2 
comes from (|13[) and cosine functions, together with the 
pole, from the expansion of the last exponent in Q - 
except of the first two j3's, which do not depend on Z 
similarly as it was done in (l5j . 

Now we should consider the AI configuration. This 
contribution may be found in an identical way and there 
is no need to repeat all steps again. It simply corresponds 
to the substitution R — > — R, so now the second term in 
brackets of {7} has the appropriate limit and contributes 
to the pole. The straightforward calculation leads to the 
similar result as above only with the modification in the 
7 matrices structure, which actually might be predicted. 



The two modes ([5]) and ([6]) exchange their roles, so the 
matrical structure is changed into (1 + 7 s ) ® (1 — 7 5 ) 
and the rest of the expression (|16p remains unaltered. 
The whole bound state pole contribution to the Green 
function is then 



iGl*+ AI (P;x,y) 



(17) 



M 2 cos(Ps/2) cos(Py/2) 2te ig 2 (0(x)+0(y » 
8tt P 2 - n 2 + ie 
x [(1 - 7 5 ) ® (1 + 7 5 ) + (1 - 7 5 ) ® (1 + 7 5 )] • 



Confronting the obtained result with that of our previ- 
ous work [15j we observe that (jTTJ) exactly cancels all 
the unwanted terms in the formula (32) and restores the 
full factorizability already on the level of the Green func- 
tion. (The absence of cos ^ in the formula (33) is only 
a typing omission.) The only apparent difference is in 9 
dependence via the two factors e 7 arising in instan- 
ton sector k = ±2 and certainly absent for the I A and 
AI configurations, which belong to k = sector. Actu- 
ally, since both factors are accompanied by (11 ±7 5 ), and 
because 7 5 (1 ± 7 5 ) = ±(1 ± 7 5 ), they simply reduce to 
e T1 ° . In consequence, appearing always in opposite pairs, 
they cancel each other. Anyway, in the massless theory, 
9 does not play any role and may be gauged away or 
simply set to 0. The above calculation reveals again the 
complexity of the field theory vacuum with underlying 
nontrivial topological structure. 
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